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ABSTRACT 

The production of 2b Al by explosions of classical novae has been computed by 
means of a hydrodynamic code that follows both the accretion and the explosion 
stages. A special emphasis has been put on the analysis of the influence of the 
initial abundances of the accreted envelope, as well as on the nuclear reaction 
rates involved. With the most recent values of chemical composition and 
reaction rates available, 26 A1 production is lowered with respect to previous 
computations. According to our results, the final contribution of novae to the 
galactic 26 A1 is at most 0.4 M 0 , which is a small part of the estimated 26 A1 in 
the Galaxy derived from COMPTEL observations of the 1809 keV emission. 

Subject headings: novae, cataclysmic variables — nuclear reactions, 
nucleosynthesis, abundances — gamma-rays :theory 
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1. Introduction 


The discovery of 26 A1 in the interstellar medium by the HEAO-3 satellite, through 
the detection of the 1809 keV y-ray line ([Mahoney et al. 1982| , [Mahoney et al. 1984|) , has 
raised the interest on potential y-ray diagnostics of nucleosynthesis in several scenarios. 
This y-ray line is produced by the decay (r ~ 10 6 years) from the ground state of 26 A1 to 
the first excited state of 26 Mg , which de-excites to its ground state level by emitting a 1809 
keV photon. This detection has been confirmed by other space missions like the SMM y-ray 
spectrometer (|3harc et al. 1985|) or several balloon-borne experiments. 


Recent measurements made with the COMPTEL instrument on-board the Compton 
Gamma-Ray Observatory (CGRO) have provided a map of the 1809 keV emission in the 
Galaxy. This map shows an extended diffuse emission along the galactic plane, with a 
peculiar large scale asymmetry about the galactic center and a clumpy structure with 
several noticeable hot spots (Pichl et al. 1995|) . The first interpretations of this map 
have suggested that novae and low-mass stars cannot be the major contributors to the 
1809 keV emission, because their low individual yields and their high frequency should 
provide a smooth distribution, in contradiction with the irregular appearance of the 
emission (see Prantzos fc Diehl 1990 for a review). The general model suggested considers 
a two-component origin of the emission: a global component following the spiral pattern of 
the Galaxy (a presumed site of massive star formation) on which several localized regions of 
intense activity, such as the Vela region, are superimposed ( [Chen, Gchrcls fc Diehl 1995|) . 
However, observations of four nearby supernova remnants with COMPTEL do not provide 
evidence for the 1809 keV emission, although the uncertainties in the distances do not allow 
to put severe constraints on 26 Al production by supernovae ([Knodlscdcr et al. 1996a|) . 
Finally, recent attempts to understand the 1809 keV map of the sky have analyzed in more 
detail its correlation with the spiral structure QKnodlscdcr et al. 1996b|) . They derive a 
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total 26 A1 mass of 2.5 M 0 , from which at least 0.7 M 0 can be attributed to massive stars, 
stressing that they cannot exclude that a large fraction of 26 A1 is produced by novae or 
low-mass AGB stars. Therefore, with this still unclear panorama, it is worth studying the 
role played by classical novae in the synthesis of galactic 26 A1 . 


As pointed out by [Ward fe Fowler 1980| , 26 A1 production requires moderate peak 
temperatures, of the order of T pea k < 2xl0 8 K, and a fast decline from maximum 
temperature, conditions that are commonly achieved in nova outbursts. In the early 80s, 
one-zone model calculations of explosive hydrogen burning nucleosynthesis by |Hillcbrandt| 
ik Thielemann 1982| and |Wicschcr ct al. 1986| suggested that classical novae might produce 
sufficient amounts of 26 Al to account for some of the observed meteoritic anomalies but 
would not represent major galactic sources. But these calculations used only solar or 
CNO-enhanced envelopes. New one-zone model calculations by [Weiss fc Truran 1990| and 
Nofar, Shaviv fc Starrficld 199 1| , on the basis of ONeMg white dwarf stars, produced large 
amounts of long-lived radioactive nuclei, such as 22 Na and 26 Al , concluding that these novae 
might represent important (though not dominant) sources of the galactic 26 A1 . However, 
according to [Nofar, Shaviv fc Starrficld 1991| , no 26 Al production results for models with 
Tpeak > 2.7x10 s K. Recent calculations have refuted this result (|Politano ct al. 1995, 


Starrficld ct al. 199B|) , demonstrating the crucial role played by convection in order to carry 
some 26 Al to the outer, cooler layers of the envelope, where its destruction through proton 
captures is prevented. But these fully hydrodynamical calculations involving ONeMg white 
dwarfs assume a composition for the white dwarf core (which is partially mixed with the 
envelope) based on old calculations of C-burning nucleosynthesis by |Arnctt fc Truran 
|1969| that need to be updated. Since the resulting nucleosynthesis is very sensitive to the 
envelope’s composition, it is important to adopt a more realistic one. On the other hand, 
there are large uncertainties on some crucial reactions involving 26 Al production. Thus, it 
is important to analyse the influence of the different prescriptions available on the final 
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26 A1 yields in the ejecta of classical novae. With this aim we have computed a series of 
hydrodynamic models of nova outbursts for several white dwarf masses and accretion rates, 
adopting updated initial chemical compositions and recent prescriptions for the crucial 
reaction rates. 


Model and results 


A one-dimensional, lagrangian, implicit hydrodynamic code has been developed, in 
order to follow both the hydrostatic accretion phase and the resulting hydrodynamic 
explosion (|Jose 1996| and [Jose fe Hernanz 1997|) . A time-dependent formalism for convective 
transport has been included ([Wood 1974]) , where partial mixing in convective shells is 
treated by means of a diffusion equation QPrialnik, Shara fc Shaviv 1979Q . With this code 
we have performed a series of computations of the complete evolution of accreting white 
dwarfs with masses ranging from 1 to 1.35 M 0 , accreting at rates between 2xl0~ 8 and 
2xlO~ 10 MQ.yr” 1 and with an initial luminosity of 1CA 2 L 0 (see |Hernanz et al. 1996 for 
an application of this code to the study of 'Li synthesis in novae). The matter transferred 
from the companion is assumed to be solar-like, but some process of mixing between 
the envelope and the outermost shells of the underlying core (shear mixing, diffusion) is 
assumed to take place. This assumption is based on the current prediction that enhanced 
CNO or ONeMg abundances are required in order to power the nova outburst and also 
to explain some of the observed abundances ( |Livio fc Truran 1994 . See also Prialnik !& 


[Kovetz 1995| and |Politano et al. 1995| for recent calculations of CO and ONeMg novae, 
respectively). The problem of the initial composition of nova envelopes is complicated and 
far from being understood. A compromise solution is to adopt a mixture with 50% by mass 
of core abundances. This may be considered as representative of the degree of mixing, 
in view of the mean metallicities observed in the ejecta of ’true’ ONeMg novae ( [Li vio fe 
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Truran 1994. See Kolb fc Politano 1996 for an analysis of the influence of the degree of 
mixing in the synthesis of 26 A1 ). The composition of the underlying core has been taken 
from recent detailed evolutionary models, specially for the case of ONeMg white dwarfs, 
which are the main contributors to 26 A1 production. These stars are made basically of 16 O 
and 20 Ne ([Dominguez. Tornambe fe Isern 1993| and [Rftossa, Garcia-Berro &; Iben 19%|) , 
since magnesium is almost absent (so we will call them ONe models). We want to stress 
that this issue has important consequences on the final yields of some species, like 26 A1 . 

A complete reaction network has been directly coupled to the code to follow the 
detailed evolution of 100 nuclear species, ranging from 1 H to 40 Ca, linked through more than 
370 updated nuclear reactions. Both the ground and the isomeric states of 26 A1 have been 
included ( 26 AD and 26 Al m , respectively). The influence of nuclear uncertainties on 26 A1 
yields has been investigated by |Coc et al. 1995| , with the use of the semi-analytical model 
of [MacDonald 1983| , for the cases considered by Politano et ah 1995| (which have a high 
initial content of 24 Mg , ~0.1 by mass). They stressed the importance of the E p = 188 keV 
resonance in 26 AD (p,y) 27 Si which is not included in the |Caughlan fc Fowler 1988| rate 
(but whose strength was measured by IVogclaar 1989| ) and also of the large uncertainties 
still affecting the 25 Al(p, 7 ) 26 Si rate. When 24 Mg is the primary source of 26 Al , formation 
of 26 AD through 24 Mg(p, 7 ) 25 Al(/?+) 25 Mg(p, 7 ) 26 AD is bypassed by 25 Al(p, 7 ) 26 Si(^+) 26 Al m , 
at high temperature. Another important reaction for 26 Al production is 2 ‘Al(p,a) 24 Mg 
which competes with 27 Al(p, 7 ) 28 Si to partially recycle the initial Mg-Al material. While 
the 2 ' Al(p, 7 ) 28 Si rate is sufficiently known in the temperature range of interest, the 
2 ' Al(p,a) 24 Mg one is still uncertain but appears to be lower than the one found in |Caughlan 
|fe Fowler 1988| . Since then, the resonances affecting this rate in the domain of temperature 
of novae have been studied by |Tinnnermann et al. 1988[ and IChampagnc et al. 1988[ Using 
two different techniques, they derive similar upper limits for the strengths which are several 
orders of magnitude lower than those used by ICaughlan fc Fowler 1988| . Instead of being 
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recycled in the Mg-Al region, the nuclear flow mainly leaks out towards higher Z elements, 
with the result of a reduced 26 A1 production. Concerning the importance of the initial 
chemical composition, we stress that the sum X( 24 Mg )+X( 25 Mg ), and not only X( 24 Mg ), 
has to be considered as a seed for 26 AH production: 24 Mg is readily transformed into 25 Mg 
by the fast 24 Mg (p, 7 ) 25 Al reaction followed by the beta decay of 25 A1, provided that the 
leak through 26 Si can be neglected. Furthermore, unlike the CNO cycle, neither the NeNa 
nor the MgAl cycles are expected to be closed and part of the initial neon and sodium 
will be transferred to the MgAl region. Indeed, in the temperature range 1.5 < T§ < 4 
the rate of the reaction 23 Na(p,ct) 20 Ne is approximately equal to that of 23 Na(p, 7 ) 24 Mg 
according to |Caughlan fe Fowler 1988[ or up to an order of magnitude lower following a 
recent reevaluation by |E1 Eicl fc Champagne 1995| . Hence, a fraction of the initial X( 23 Na) 
will reach the MgAl region providing additional seeds while the rest ends up as 20 Ne. A 
fraction of 22 Ne can also be transformed into 23 Na. At the temperatures prevailing in nova 
outbursts, the rate of 20 Ne(p, 7 ) is low and most of the initial 20 Ne remains unburned. 
However, in classical novae involving massive white dwarfs, a significant fraction of 20 Ne 
can be transformed through the NeNa cycle and via 23 Na(p, 7 ) 24 Mg into 24 Mg. Accordingly, 
even in the absence of magnesium in the initial composition, 26 A1 can be produced from 
23 Na, 22 Ne or even some 20 Ne in the most extreme cases. 


Snapshots of the evolution of several isotopes relevant to the 26 A1 synthesis ( 23 Na , 
24 Mg , 25 Mg , 26 Al 3 , 26 Al m and 27 A1 ) are shown in figure 1, for a 1.25M Q ONe nova, with 
an accretion rate of 2xlO~ 10 M 0 .yr _1 . The upper panel corresponds to the time when 
the temperature at the burning shell has already reached 10 8 K. Nearly flat abundance 
profiles along the accreted envelope are found, since the convective timescale is similar to 
the characteristic timescale of most of the dominant nuclear reactions. Proton capture 
reactions on 24 Mg and 25 Mg provide the main source of the Al-group nuclei which increase 
their initial abundances. When the temperature in the burning shell rises up to 2.1x10 s K 








(second panel), the nuclear timescale becomes much shorter than the convective timescale. 
Hence, only partial mixing between adjacent convective shells occurs and non-flat profiles 
are obtained. At this time, the abundance of 26 A1 in the burning shell attains its maximum 
value, X( 26 A1 ) ~ 7.7xl0 -3 by mass. 26 A1 is being progressively destroyed near the burning 
shell (panel 3) due to proton captures, since the seed nuclei 23 Na , 24 Mg and 25 Mg have 
almost been exhausted. The final abundances in the ejecta are shown in the lower panel, 
which corresponds to the time when the white dwarf envelope has already expanded to 
a size of R w d ~ 10 12 cm. The amount of 26 A1 in the envelope reflects the efficiency of 
convective transport, which carries 26 A1 isotopes to the outer, cooler layers of the envelope, 
where they can avoid destruction through (p,y) reactions. 


The 26 A1 yields obtained in all models computed are summarized in table 1. All 
26 A1 production appears in its ground state, since the isomeric state doesn’t play any 
role. ONe novae are more important producers than CO ones, due to the fact that peak 
temperatures attained during CO nova outbursts are not high enough to break the CNO 
cycle. Therefore, the production of 26 A1 comes from the initial 23 Na , 24 Mg and 25 Mg . It is 
worth noticing that the initial abundance of 24 Mg is much lower than the one adopted in the 
calculations performed by [Tolitano ct al. 1995 , which were based on the old nucleosynthesis 
computations of carbon-burning by [Arnett &; Truran 1969| . However, according to the 
above discussion, the presence of 25 Mg , 23 Na and 22 Ne partially compensates this absence. 


Only some combinations of peak temperatures around T peak ~ 1-2 xlO 8 K and rapid 
evolution from maximum favor 26 A1 generation. In this sense, explosion in low-mass 
ONe white dwarfs, where the lower gravity provides softer explosions with lower peak 
temperatures, favor 26 Al production. Our calculations (see table 1) show a minimum in the 
26 Al production near M wc j = 1.25 M 0 , which fairly agrees with the trend obtained in the 
semi-analytical analysis of Coc et al. 1995| , but doesn’t follow the monotonic decrease in 
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the 26 A1 production as the white dwarf mass increases reported by Politano et ah 1995 . We 
point out, however, that the amount of 26 A1 ejected into the interstellar medium decreases 
as the white dwarf mass or the mass accretion rate increase, due to the fact that the mass 
of the envelope is smaller. It is also worth noticing that the amount of 26 A1 present in the 
ejecta of the models described in this paper is systematically lower than the mean values 
obtained by [Politano ct al. 1995| . Two main effects influence this result: first, a different 
initial composition (mainly the lower initial content of Mg in our calculations), and second, 
the specific prescriptions adopted for the reaction rates of interest for the 26 Al synthesis. 
The effect of the nuclear reaction rates has been checked by means of 2 models of 1.25 M 0 
ONeMg white dwarfs, accreting mass at a rate 2xlCP 9 M 0 .yr -1 , which are similar to 
one of the models analyzed by Politano et al. (1995). Both models differ in the specific 
rates adopted for the reactions 26 A1 (p,y) 27 Si and 2 'Al(p,a) 24 Mg. Since these reactions 
play no role during the accretion phase, the total amount of mass accreted as well as the 
proper pressure at the onset of the explosive phase are the same in both models. Hence, 
the differences obtained in the final yields reflect exactly the different prescriptions for the 
reaction rates adopted. When the rates from Vogelaar (1989) and Champagne et al. (1988) 
are used, we obtain a mean value of X( 2b Al ) = 1.4 xl0~ 3 by mass in the ejected shells, 
whereas X( 26 A1 ) = 4.8 xlO” 3 is found when the Caughlan and Fowler rates (1988) are 
used (|Politano et al. 1995| obtained X( 26 A1 ) = 9.45 xlCF 3 by mass, for a similar model). 
In addition, the ratio 26 A1 / 2 'A1 increases from 0.4 to 3.0 in the second model. This result 
is a direct consequence of the differences between the two prescriptions for the rate of 
2 ‘Al(p,a) 24 Mg: the lower cross-section reported by Champagne et al. (1988) results in a 
less efficient recycling of 2 'A1 into 24 Mg , which essentially will end up in the form of 26 Al . 
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3. Discussion and conclusions 


The production of 26 A1 by classical novae is very sensitive to the initial composition 
of the envelope and to the nuclear reaction rates adopted. ONe novae are more important 
26 A1 producers than CO ones, because seed nuclei for the NeNa and MgAl cycles are almost 
absent in the latter ones. For the same reason, the amount of 26 A1 synthesized in ONe 
novae depends on the initial composition of the white dwarf core. When using the recently 
available chemical abundance profiles by |Kitossa, Garcia-Berro fc Iben 1 0T7B| , a lower 26 A1 
production is obtained. Furthermore, some improvements in the nuclear reaction rates since 
Caughlan fe howler 1988| also lead to a lower 26 A1 production. 


The amount of 26 A1 injected into the interstellar medium by ONe nova events decreases 
as the mass of the underlying white dwarf increases. Hence, low-mass white dwarfs are 
most likely candidates for 26 A1 production. Contribution from low-mass white dwarfs is 
favored both by the higher 26 A1 production and the higher ejected mass. But white dwarfs 
of masses lower than ~1.1 M 0 are expected to be CO white dwarfs, which are unable to 
produce important quantities of 26 A1 (see table 1). Also, the observations of some Ne-novae, 
such as QU Vulpecula 1984, indicate high ejected masses (~ 10” 3 M 0 ), unobtainable if a 
massive ONe nova is the responsible of the explosion. This has led some authors to propose 
a new scenario (EE ara 1994] , |Shara fc Priainik 1994| ): low-mass CO white dwarfs, undergoing 
episodic accretion phases at high rates, experience large metal enrichments from the ashes 
of He-burning. The subsequent nova explosion, correlated with a phase of a lower accretion 
rate, produces at the same time high 26 A1 abundances and high ejected masses. But these 
models require a fine tuning of some parameters relative to the cataclysmic variable, such 
as the mass accretion rate. 


A crude estimate of the total production of 26 A1 by novae can be computed by means 


of the expression from |Wciss fc Truran 1990|: 
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M / 26 An _ n A X( 26 A1) frac(Ne) M ej (M 0 ) R nova (events yr *) 
^ ’ ' 4 2x10" 3 0.25 2xl0 -5 40 


where X( 2b Al ) represents the mean mass fraction of 26 A1 in the ejecta, frac(Ne) the fraction 
of Ne-novae over the total number of classical nova outbursts (between 0.25 and 0.53, 
but typically ~ 1/3, see [lruran fc Livio 1986| , [Livio fc Truran 1994]) , M e j the amount of 
mass ejected in an outburst, and R llova the nova rate (~ 46yr _1 , |Higdon fc lowlcr 1987| 
or ~ 20yr _1 , Pella Valle fc Livio 1994|) . Adopting these estimates of R nov a and frac(Ne) 
and our most favorable ONe nova case, for which M e j ( 26 A1 ) = 1.7 xlO -8 M 0 , we derive 
a maximum contribution of classical nova outbursts to the galactic 26 A1 in the range 0.1 
to 0.4 M 0 . Although higher abundances of 26 A1 can be obtained for models with higher 
enrichment from core material (i.e., 75%), these may be considered as very extreme cases, 
in view of the typical metallicities observed in the ejecta of classical nova outbursts. In 
summary, the contribution of novae to galactic 26 A1 is small as compared with the one 
required to explain COMPTEL measurements (between 1 and 3 M 0 ). This is in agreement 
with the most accepted hypothesis of young progenitors as sources of galactic 26 A1 . 


We want to stress that the two aspects that are crucial for the final 26 A1 yields are 
far from being understood; these are the process of mixing between the core and the 
envelope, which determines the initial chemical profile, and the exact rates of the nuclear 
reactions involved in 26 A1 synthesis. This latter topic will surely be improved with the new 
compilation of nuclear reaction rates ([NACRE 19971 ). 


We are specially indebted to J. Isern, for valuable discussions about the topic of 
this paper. We acknowledge for partial support the CICYT (ESP95-0091), the DGICYT 
(PB94-0827-C02-02), the CIRIT (GRQ94-8001), the AIHF 95-335, the Human Capital and 
Mobility Programme (CESCA/CEPBA) and the PICS 319. 
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Table 1: 26 A1 yields and ejected masses for some nova models 


Comp. 

M wd (M 0 ) 

M^yr- 1 ) 

X( 26 A1 ) 

M.,( 26 A! )(M S ) 

CO 

1.0 

2xlO -10 

1 .6xl0~ 5 

3.5xlO -10 

CO 

1.15 

2xlO~ 10 

4.7xl0 -5 

6 .2xlO -10 

ONe 

1.15 

2 xlO~ 10 

9.3xl0" 4 

1.7xlCT 8 

ONe 

1.25 

2 xlO~ 10 

5.4xl0” 4 

7.7xlCT 9 

ONe 

1.25 

2 xl0~ 8 

5.7xlCT 4 

4.7xlCT 9 

ONe 

1.35 

2xlO~ 10 

7.2xl0 -4 

3.1xl0 -9 
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Fig. 1.— Snapshots of the evolution of 23 Na , 24 Mg , 25 Mg , 26 Al 9,m and 2 ' A1 along the 
accreted envelope, for a 1.25 M 0 ONe novae accreting at a rate M = 2 xlO~ 10 Mg.yr^ 1 . 
The successive panels, from top to bottom, correspond to the time for which the temperature 
at the burning shell reaches 10 8 , 2.1 xlO 8 , T max = 2.44 xlO 8 , plus an additional case for 
which the white dwarf envelope has already expanded to a size of R w( j ~ 10 12 cm. The arrow 
indicates the base of the ejected shells. 






